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PREFACE 


A previous report. Reference 1, prepared under this contract, described 
a computer program, designated RETSCP, for the analysis of rocket 
engine thrust chambers with cyclic plasticity. That report included a 
complete description of the program logic, together with a users manual. 
It is the purpose of this report to illustrate the detailed application of 
the RETSCP program. Volume I describes the analysis of a re- 
generatively cooled OFHC copper combustion chamber. Operating 
conditions are such that plastic strains dominate and fatigue life 
is in the low cycle regime. Volume II describes the analysis of an 
attitude control thrust chamber. Operating conditions for that 
engine are such that the material behavior is elastic, or on the 
threshold of the elasto-plastic regime; thus, the fatigue life is in 


the high cycle regime 


SUMMARY 


A three-dimensional finite element elasto-plastic strain analysis was 
performed for the throat section of a regeneratively cooled rocket 
combustion chamber. The analysis employed the RETSCP finite 
element computer program (NASA CR-1 34640). The analysis Included 
thermal and pressure loads , and the effects of temperature dependent 
material properties, to determine the strain range corresponding to 
the chamber operating cycle. The analysis was performed for chamber 
configuration and operating conditions corresponding to a hydrogen- 
oxygen combustion chamber which was fatigue tested to failure at 
the NASA Lewis Research Center. The computed strain range at 
typical chamber operating conditions was used in conjunction with 
OFHC copper isothermal fatigue test data to predict chamber low-cycle 


fatigue life. 


1,-s ^ 


INTRODUCTION 

A new generation of htgh performance liquid rocket engines is being 
considered for main propulsion in Space Transportation Systems . The 
high performance goal which has been established for these engines 
demands high chamber pressures with resulting high heat flux levels 
in the combustion chamber. For Space Transportation System applications, 
engine reusability is a prime objective. Thus, the stress analyst must be 
able to define the life potential of a given design under conditions of 
thermal and pressure cycling. 

Combustion chambers can be designed to operate at suitable temperature 
levels by emoloying regenerative cooling . The combination of high heat 
flux and high pressure necessitates the use of materials with high thermal 
conductivity. Such materials permit practical values of wail thickness 
without creating excessive hot gas side wall temperatures. Typical 
construction consists of a copper based alloy combustion liner, containing 
milled coolant passages, combined with a nickel alloy structural jacket. 

The thermal analysis of combustion chambers such as this follows established 
methodology. Empirical data is available for estimating combustion gas 
and coolant heat transfer. The temperature distribution within the structure 
is computed using one of many available finite difference computer codes 
for transient or steady state heat conduction in two or three dimensions. 

The state of stress in regeneratively cooled rocket chambers varies in 
three dimensions. The situation is further complicated by the strong 
dependence of material properties on temperature. It is the range of strain. 
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through which a given chamber is cycled, that determines fatigue life. 

A numerical method of structural analysis must be employed for such 
complex configurations . 

The objective of Volume I is twofold. First, there is the application 
of a three dimensional finite element computer program to determine the 
thermal strain behavior of a specific regeneratively cooled copper 
combustion chamber. An indepth study of this type, not only provides 
strain range data for fatigue life analysis; but also lends insight 
into the structural response of the chamber under complex loading 
conditions . The copper chamber configuration and operating conditions 
correspond to a chamber which was fatigue tested to failure at the 
NASA Lewis Research Center. 

The second objective is to illustrate in detail the use of the RETSCP 
finite element computer program for chamber configurations where 
plastic strains dominate; that is, low cycle fatigue. Thus, considerable 
attention is given to the finite element model, material properties, and 
execution of the program. 

The combustion chamber configuration, operating conditions, and test 
results are described in sufficient detail to conduct the strain analysis 
and compare results. The strain analysis is then presented in such 
a manner that the work could be extended to other configurations and 
operating conditions. The fatigue life analysis section describes the 





predictive method used for the copper chamber; and, also indicates 
methods which could be employed for other studies . The 
Concluding Remarks include discussion of future applications for 
the RETSCP computer program, including refinements for fatigue 
life prediction, and phenominalogical studies of the thermal 
strain behavior. 
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COMBUSTION CHAMBER CONFIGURATION 


The fatigue life analysis was conducted for the throat section of a 
specific regeneratively cooled rocket thrust chamber which was 
fatigue tested to failure at the NASA Lewis Research Center. Since 
the experimental results are unpublished, those facets of the 
geometry, operating conditions, and test data which pertain to the 
life analysis are given below. 

Chamber Geometry 

The overall combustion chamber geometry and construction data are 
shown in Figure 1. The basic component of the chamber for high 
heat flux and high pressure is the OFHC copper liner. This liner 
contains sixty milled axial coolant flow passages. The high thermal 
conductivity of the copper produces two beneficial effects. First, 
from a manufacturing standpoint, the material permits practical 
values of wall thickness. The allowable hot wall operating temperature 
together with the imposed heat flux value, determine the wall thickness 
Thus, for materials with low thermal conductivity, such as nickel 
alloys, extremely thin walls are required. The second effect, due to 
high thermal conductivity, is that the ribs which divide the cooling 
passages are very effective as cooling fins. 



STATION 



Figure 1 . Combustion chamber geometry 



Chamber structural integrity is maintained by the electroformed 
nickel closeout jacket. The nickel has low thermal conductivity 
(as compared with copper); and, therefore, low temperatures are 
encountered throughout the operating cycle . Thus , the nickel 
jacket restrains the chamber pressure hoop and thrust loads during 
engine operation. Note, that the nickel structure contains integral 
coolant manifolds; and , that the coolant flow direction is from 
chamber nozzle toward engine injector. 

The chamber is instrumented to measure chamber and coolant pressures 
as well as coolant inlet and outlet temperature. At the throat section, 
the nickel jacket and copper liner temperatures are measured at four 
circumferential stations. Nickel jacket thermocouples are attached 
to the outer wall surface. Copper liner thermocouples are located 
nominally at 0.050 in. (0.127 cm) from the hot gas side surface, 
circumferentially centered within the liner ribs . One rib temperature 
thermocouple was inoperable due to damage during chamber fabrication. 

Operating Conditions 

The fatigue test thrust chamber is a regeneratively cooled hydrogen- 

oxygen chamber which is cooled by cold hydrogen near 50°R (28°K) 

6 , 2 

inlet temperature and 1066 psia (7.35x10 N/m ) inlet pressure. 
Coolant flow preceeds chamber firing and by-passes the chamber; 
thus, producing precool to a uniform temperature near 50°R (28° K). 




The thrust chamber pressure pulse Is Indicated by the solid line curve 

in Figure 2. During the chamber start transient, the chamber pressure 

6 2 

is increased to 610 psia (4.21 x 10 N/m)ina2.25 second interval. 

The step transient start is required by facility limits. Measured 
thrust force durim the 2.0 second operating phase was 4686 lbs. (20840 N). 
Thrust chamber shutdown occurs over a 1.75 st jond interval. The 
hydrogen coolant bypass flow is maintained for the 2.1 second period 
between pulses. The test chamber was subjected to pulse trains 
consisting of one to eleven pulses of the type described. 

6 2 

Combustion gas conditions at the throat station are 347 psia (2.39X10 N/m ) 

static pressure and 6143°R (3413 °k) adiabatic wall temperature. Coolant 

temperature and pressure at the throat during steady state firing are 

6 2 

132°R (73°K) and 1045 psia (7.21 x 10 N/m ) respectively. 
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Test Results 

The measured temperature response at the chamber throat is depicted 
by the dashed line curve ( A T) in Figure 2 . The value A T represents 
the difference between the measured copper liner temperature and the 
nickel jacket outer wall temperature. During the engine start transient, 
the copper liner heats at a greater rate than the low conductivity 
nickel jacket. After the copper has reached steady state temperature, 
the nickel approaches steady state conditions; thus , generating the 
thermal overshoot in Figure 2 . A similar cooling overshoot occurs 
during thrust chamber shutdown. 

The thrust chamber was tested at the operating conditions previously 
described for pulse trains of 1,1,1,1,5,4,3,11,1,2,3,5,3,2 pulses 
each for a total of forty-three cycles . Visual inspection afte- 
forty-third pulse revealed a thru-crack from the coolant passa. '.w the 
combustion chamber interior at the chamber throat . The thru-crack 
probably occured on the thirty-ninth cycle when a change in the test 
cell noise pattern was perceived. 

It should be noted that geometric changes were observed along the 
combustion chamber length as the test program progressed, coolant 
channel bulging and surface roughening of the hot gas side copper 
surface occurred. These effects were apparent after only a few cycles 
of testing. Associated with these changes, the thermal characteristics 
of the engine were altered, and the measured copper temperature levels 
showed a gradual increase with each successive firing. 




STRAIN ANALYSIS 

It is well known that fatigue life depends on the cyclic strain range 
magnitude. One could compute the strain range for each loading 
cycle using the RETSCP program. Reference 1, and determine the 
amount of fatigue damage by applying the linear cummulative damage 
law, Reference 2. In th interest of computional practicality, the 
strain analysis is based on a set of typical or average material properties. 

Finite Element Model 

The RETSCP strain analysis employed the thirty-four element model 
of the copper chamber throat section shown in Figure 3. The cross 
section shown is one of 120 similar symmetry sections. Each element 
has the same thickness of 0.010 in. (0.025 cm) in the z - direction. 

Of course, the throat station was chosen for analysis since it is the 
location of the extreme thermal condition. 

The number of elements was selected to adequately represent the 
structure, More elements were used in locations where large changes 
in strain and temperature are expected; for example, near the hot gas 
side surface. It was noted in Reference 1, that the RETSCP program 
is based on an isoparametric element; and thus, few elements are 
required to model the structure. 



(0., .340) (.0859 .3346) 



(a) Geometry and Element 
Numbers 


(b) Node Numbers and 
Partition Numbers 


Figure 3. Finite element model. 
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Boundary conditions in the RETSCP program are in the form of prescribed 

nodal point forces or nodal point displacements. The symmetry of the 

34 - element model was previously noted. The corresponding boundary 

condition is zero displacement of boundary points normal to the symmetry 

plane, and free movement of nodal points in the radial direction. This 

condition is indicated by the roller notation in Figure 4. The symmetry 

o 

condition along the 3 surface is accommodated by a coordinate 
transformation within the RETSCP program. 

The surfaces exposed to the coolant and combustion gasses experience 
prescribed force in the direction normal to the surface. The pressure 
is replaced by force components at each nodal point as indicated in 
Figure 4. Each point carries load from adjacent pressure surfaces. 

Of course, the force components on the combustion surface are zero 
for engine-off conditions . 

The boundary conditions in the z - direction require special consideration. 
A very thin section is taken at the throat, so that the rocket engine 
contour curvature can be neglected. That is, constant thickness elements 
are used. Under load, the section remains plane; and also, the total 
load on the section is equal to the engine thrust load. This corresponds 
to the generalized plane strain condition. 


Symmetry plane 


Symmetry plane 



Figure 4. Force and displacement (symmetry) boundary 
conditions . 
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Recall that in the RETSCP program, either forces or displacements are 
prescribed at the nodal points. In this case, we prescribe zero 
displacement of all points on one face , and a single displacement 
value for all points on the opposite face . Since the nodal point forces 
are output by the RETSCP program, we make a posteriori check that 
the sum of the forces in the z - direction equals the thrust force. 

If not, the z - displacement is adjusted and the process repeated. 



The initial displacement estimate is based on the generalized plane 
strain equation applied to the nickel closeout structure; namely , 


Reference 4 , (Symbols are defined in Appendix A) 


1 

e z = T 


°z- v ( °x + 0 y) 


( 1 ) 


Stresses in the y- direction are neglected. The x - stress in the 
nickel is approximated by the pressure hoop stress. That value is 
adjusted due to the thermal expansion force exerted by the copper liner 
in the hoop direction. The copper is at yield stress conditions. 
Similarly, the z - stress is approximated by the axial thrust stress 
which is corrected due to the force exerted by the copper in the axial 
direction at yield stress conditions. The computed strain is then 
converted to displacement, and adjusted by the average thermal growth 
of the nickel to give the initial prescribed displacement estimate. 


* 




i 

i 
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Thus, the terms in equation (1) are: 


internal 

pressure 


throat 

diameter 


average copper 
. yield stress . 


nickel radial thickness 


copper radial thickness 
at coolant passage 


+ average copper copper cross 
thrust yield stress sectional area 

' . m i m, i mm — — " — ■ » 

[nickel cross sectional areal 


In principle, the above procedure will lead to a consistent set of results 
with the proper axial force solution. In the case of the copper chamber, 
the copper strain behavior is so dominated by the thermal growth , that 

the small axial displacement value has little effect. Thus, correction of 

* 

the initial guess value is not warranted; and, although the stress field in 
the nickel may be somewhat erronious , the strain range of the copper 
liner is not effected. 


Temperature Distribution 


A complete thermal analysis was performed for the regenera tively cooled 
copper combustion chamber. The calculation of heat transfer coefficients , 
overall heat balance , and coolant properties was performed by NASA and 
is furnished data. 


Transient heat conduction temperature distributions were computed for an 
89 node model of the configuration shown in Figure 3 , by using a 
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finite difference computer program. Based on the computer output, 
interpolation was used to obtain the average temperature for each 
of the 34 finite elements , Figure 3 . 

The predicted temperature values, corresponding to the measurement 
locations , showed excellent agreement with the values measured 
during the first few cycles of testing. As previously noted, the 
measured temperature values showed a gradual increase as testing 
progressed. Thus, the temperature distributions were recomputed 
with input boundary conditions such that results would agree with 
temperature values from the twentieth test cycle. The twentieth cycle 
is approximately one half of the chamber life; and, is thus assumed 
to represent average or typical conditions . Input conditions were 
selected in the program such that computed results agreed with the 
thermocouples in the copper liner. Measured values on the nickel 
were somewhat lower than the recomputed values which indicated that 
convection to the surroundings and axial conduction heat transfer 
were significant. The nickel elements were further adjusted to 


account for this deviation 




The NASA furnished sets of temperature data for the finite element model 
are listed in Table I, for the four critical points of the operating 
cycle indicated in Figure 2; that is; (a) precool-intercool, (b) start 
transient peak overshoot, (c) steady state chamber firing, and 
(d) cooldown transient peak overshoot. Sample results, taken from 
the four sets of 34 temperature difference values are shown in Figure 5. 
Note that the maximum temperature value is 1493°R (829° K) which 
corresponds to the temperature difference of 993°R (552° K), which 
occurs during the steady state chamber firing. Temperature differences 
in Table I and Figure 5 are in units ( °R) suitable for the RETSCP program. 
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TABLE I 

* 

NASA FURNISHED TEMPERATURE DIFFERENCE (A T = T - T ref) 
FOR TWENTIETH TEST CYCLE (T ref - 500°R). 


Element 

A Ta 

(t=0 or 8.1 sec) 

A Tb 

(t= 2.4 sec) 

A Tc 

(t= 4.2 sec) 

A Td 

(t= 6.0 sec) 

1 

-450. 

-153 

+92 

-135 

2 

-450. 

-153 

+95 

-135 

3 

-450. 

-153 

+95 

-135 

4 

-450. 

-30 

+127 

-219 

5 

-450. 

+18 

+170 

-224 

6 

-450. 

+37 

+189 

-224 

7 

-450. 

+235 

+330 

-264 

8 

-450. 

+249 

+344 

-262 

9 

-450. 

+316 

+401 

-278 

10 

-450. 

+327 

+413 

-275 

11 

-450. 

+483 

+557 

-290 

12 

-450. 

+493 

+568 

-288 

13 

-450. 

+601 

+671 

-295 

14 

-450. 

+606 

+677 

-293 

15 

-450. 

+830 

+893 

-303 

16 

-450. 

+814 

+877 

-302 

17 

-450. 

+776 

+841 

-300 

18 

-450. 

+727 

+794 

-297 

19 

-450. 

+708 

+776 

-295 

20 

-450. 

+847 

+905 

-302 

21 

-450. 

+839 

+903 

-302 

22 

-450. 

+799 

+864 

-300 

23 

-450. 

+753 

+820 

-297 

24 

-450. 

+743 

+812 

-296 

25 

-450. 

+885 

+947 

-302 

26 

-450. 

+878 

+941 

-302 

27 

-450. 

+846 

+911 

-300 

28 

-450. 

+817 

+882 

-298 

29 

-450. 

+809 

+875 

-298 

30 

-450. 

+931 

+993 

-304 

31 

-450. 

+912 

+975 

-302 

32 

-450. 

+896 

+9 59 

-301 

33 

-450. 

+876 

+941 

-300 

34 

-450. 

+864 

+930 

-300 


♦Note: Values are in °R (1.8 x°K) which is suitable unit for RETSCP. 
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Figure 5. Computed temperature difference 
(AT»T-T ref) for twentieth test cycle (T ref = 500°R) . 

Values are in °R (1.8 x °K) which is suitable unit for RETSCP. 
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Material, Properties 

Materials properties are input as functions of temperature in the RETSCP 
program. Also, RETSCP plasticity analysis is based on bi-linear 
representation of the uniaxial test stress strain curve. It is well known 
that materials such as OFHC copper, work harden, or work soften, or 
harden then soften upon plastic strain cycling, Reference 5. 

Experimentally determined stress-strain curves are seldom available 
as functions of temperature and the number of strain cycles. Thus, it 
is necessary to construct stress-strain curves based on theory together 
with the best available experimental data. 

It turns out that for plastic behavior, the true stress-true strain curve 
appears as a nearly straight line when plotted on logarithmic coordinates , 
Reference 6. Room temperature data for annealed OFHC copper, Reference 7, 
is presented by the open circles in Figure 6. Other curves were constructed 
end are shown in Figure 6, based on the 0.2 % offset yield stress, and the 
ultimate tensile stress at 40% strain which is assumed to correspond to 
the start cf necking for annealed OFHC copper as reported in Reference 8. 
These data all apply for virgin material. 

The cyclic stress-strain curve, Reference 9, has great value for strain 
cycle calculations . Such curves are constructed with the aid of fatigue 
test data. The test strain range, together with the measured stress range 
which occurred at one-half the number of cycles to failure, defines a 
point on the cyclic stress-strain curve. 



Figure 6. True stiess - true strain curves for OFHC copper 






The cyclic stress-strain curve for annealed OFHC copper, shown in Figure 6, 
is based on copper fatigue test data, Reference 10, at 1000°F (538°C) 

(closed circles). The effect of cyclic work hardening is quite clear 
(about 30%). Since strain hardening of copper occurs in the first few 
cycles, work hardening is accounted for by using the cyclic stress strain 
curves . 

The stress-strain data are replotted in Figure 7 for strain values up to 3%. 
The plastic portion of the curves were constructed directly from Figure 6 . 

The pre-yield portion of the curves, in Figure 7, are faired curves to 
smoothly connect the linear elastic curve, based on modulus, to the 
plasticity data. The cyclic data in Figure 7 shows 30% hardening 
(stress increase over virgin material) at 3% strain. 

Cyclic stress strain curves at other temperatures were constructed by 
increasing the stress 30% at 3% strain and redrawing the curve to the 
same shape. Results are shown by the dash line curves in Figure 8. 

The bi-linear representations, which were used for the rocket chamber 
strain calculation, are shown by the solid lines in Figure 8. Stress-strain 
data at intermediate temperature values is computed within the program 
by linear interpolation between the values shown in Figure 8. 

Note, that for small elasto-plastic strain values (less than 0.5%), the 

bi-lincar approximations in Figure 8 do not approximate the actual curve 

very well. Thus, a set of bi-linear curves for strains up to 0.5% were 
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Yield stress at -320°F (-195 C), Reference 11 
Yield stress at room temperature , Reference 10 

Yield stress at 1000°F (538^), Reference 10 


Stress at Nf/2 for strain shown, at 
1000°F (538^), Reference 10 



Figure 7 . Stress-strain curves for OFHC copper 
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10 , 000 + 

(6895) 
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Percent 

Strain 


Bi-linear cyclic stress-strain curves for OFHC copper 
up to 3 percent strain . 






constructed, and are presented in Figure 9. Si ecit. ally, these curves 


were used in the analysis for element number 7,8,9 and 10. The 
applicability of this strain level for those elements will be apparent later. 

A summary of all data used in the RETSCP program for the OFHC copper combustion 
chamber is given in Table II. 

Structural Analysis 

First, let us qualitatively examine the structural behavior of the combustion 
chamber as it is subjected to operating conditions and thermal loading 
such as depicted in Figure 2 and 5. 

o o 

The chamber is initially subjected to a 50 R (28 K) precool temperature in 
combination with the coolant pressure load. Since the copper has a larger 
coefficient of thermal expansion (contraction) than the nickel, precool 
tensile stresses result in the copper and precool compressive stresses 
are induced in the nickel jacket. The high strength nickel jacket is 
sufficiently massive to remain elastic; however, the copper stress level 
may be high enough that yielding of the copper occurs . 

Upon heating to the start transient peak overshoot level, the situation 
is reversed . The large amount of thermal growth of the copper liner 
generates plastic compressive deformation. The nickel structure behaves 
elastically in tension to resist the chamber pressure hoop loading. 
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TABLE II 


SUMMARY OF MATERIALS PROPERTY DATA 


Property 

English 

Units 

(RETSCP) 

Metric 

Units 

Thermal Expansion Coefficient: 

Nickel 

7.2xlO~ 6 in/in-°F 

13.0xl0"®m/m-°C 

Copper 

9.8xlO' 6 in/in-°F 

17.6xl0“ b m/m-°C 

Modulus of Elasticity: 

Nickel 

30.0x10^ psi 
16.6x10*1 psi 

20.7xl0*J°N/m* 
11.4x10 lu N/mf 

Copper-cold 

Copper-hot 

lO.OxlO* 6 psi 

6.9xlO +10 N/ni 

! 

Poisson’s Ratio: 



Nickel 

.31 

.31 

Copper 

.33 

.33 

Yield Stress: (p y * o G -X» T) 

Nickel 

55,000 psi 

+ 6 . 2 
379. xlO N/m 

Nickel 

0. 

+6 2 

Copper Large Strain ( °o) -hot 

9,200 psi 

63.4x10 N/m 

5.75xl0 +4 N/m^-°C 
4.14xl0 +4 N/m“ 
3.26xl0 +4 N/mJ-°C 
11.7xl0 +7 N/m Z 

Copper Large Strain ( Xi) -hot 

4.63 psi/°F 

Copper Small Strain (oo) -hot 

6,000 psi 

Copper Small Strain ( Xi ) -hot 

2.63 psi/°F 

Copper Large Strain ( oo ) -cold 

17,700 psi 

Copper Large Strain ( Xi ) -cold 

0. 


Copper Small Strain ( oo ) -cold 

11 ,000 psi 

75.8x10 0 N/m 4 

Copper Small Strain ( Xi ) -cold 

0. 

0. 

-3 -6. 

Plastic Modulus Ratio: (m=m o xl0 -X*TxlO ) 

i 


Copper Large Strain ( mo) -hot 

41.7 

33.1 V 0 F 

41 .7 . 

59.6 V °C 

Copper Large Strain ( X* ) -hot 

Copper Small Strain ( mo) -hot 

117.7 V°F 

211.9 V°C 

Copper Small Strain ( Xj ) -hot 

6.32 

6.32 

Copper Large Strain (mo ) -cold 

31.7 

31.7 

Copper Large Strain ( X 2 ) -cold 

0. 

0. 

Copper Small Strain (mo) -cold 

132. 

132. 

Copper Small Strain ( X 2 ) -cold 

0. 

0. 

Nickel (mo) 

1.0 

1.0 

Nickel ( Xa ) 

0 . 

0 . 
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The stress-strain condition at steady state chamber firing may introduce a 
variety of effects . Just as the nickel-to-copp*- r temperature difference Is 
relaxed; so is the thermal strain level relaxed. The important poir.t N that 
this relaxation takes place along the elastic line. This point will be 
illustrated later. The strain relaxation may terminate either with t ns lie 
or compressive elastic stresses in the copper; or, relaxation may be 
sufficiently great to generate tensile plastic deformation. 

The cooldown transient peak overshoot corresponds to the maximum strain 
state under cold conditions . As the chamber returns to uniform intercool 
conditions, 50°R (28°K) with coolant pressure loading, relaxation type 
behavior will again be exhibited . 

The second test cycle will produce similar qualitative behavior with two 
major differences. First, the residual strain values, for all segments of 
the cycle , will differ from the first loading cycle . Second , the material 
will have work hardened as a result of deforming plastically. A few loading 
cycles will "shake-down" the effect of the residual strain starting point. 
Also , strain hardening of the copper will diminish after a few cycles , 
Reference 12 and 13. For the present calculation, work hardening is not 
accounted for on a cycle by cycle basis; rather, hardening is included 
through the use of cyclic stress strain data. Thus, a repeatable hysteresis 
strain cycle will eventually be established. It is noteworthy, that, the 
maximum strain range for any cycle will be determined by the peak to peak 
variations in nickei-to-copper temperature difference . 



RETSCP Program Execution 

The detailed copper chamber strain cycle was computed using the 34 
element model shown in Figure 3. In the absence of cycle uy cycle 
strain hardening data, the cyclic stress-strain curves presented in 
Figures 8 and 9 , and discussed previously, were employed. Otht. 
input material properties are listed in Table II. Pressure iouds ad 
boundary conditions correspond to Figure 4. Thermal loads were 
based on the twentieth cycle temperatur e differences given in Table I 
and Figure 5 . Thus , the calculated strain cycle corresponds to an 
average or typical test cycle. 

Each computer run consists of strain analysis based on one set of 
operating conditions together with one set of prestrain data. The 
calculation starts by introducing the precool condition (a) from Table I 
together with cold material properties data and zero chamber pressure. 
The RETSCP program output consists of a stress-strain listing plus 
punch card residual strain data for each element. The punch cards are 
entered with condition (b) temperature in Table I; and, using the hot 
material properties data and maximum chamber pressure values, the 
heating transient part of the strain cycle is computed. Again, punch 
card residual strain data is output for subsequent calculations . 

A typical set of input data for the chamber strain analysis is given in 
Appendix B. The particular data is for the transient heating condition 
(b), see Reference 1 for data card content and format. 
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Input for other conditions is similar to that listed in Appendix B 


For the RETSCP program, a secant modulus plasticity iteration, 

Reference 1 , is performed within each run . Five such iterations 
were required for convergence of all elements with the bi-linear 
stress-strain curves. Total computing time for each run was 27 
minutes on the IBM 7094 computer. 

Results 

A typical computed effective stress -effective strain diagram is shown 
in Figure 10 for a specific element; namely element 34, which is the 
element at the hot gas side surface benea*n the rib, with temperature 
differences designated (864, 930, -300, -450) in Figure 5. Figure 10 
depicts six strain loop calculations . The aforementioned effect of 
different pre-strain values on successive loop results is apparent. The 
net effect of strain hardening is automatically included by the use of 
cyclic-stress strain data. The important point is that a stable hysteresis 
loop is achieved which is representative of the strain range for that 
particular element on the twentieth cycle of thrust chember operation. 

Referring to Figure 10, consider the detailed behavior of that particular 
element. Initial pre-cooling of the virgin material produces tensile 
yielding to the point (a). During the start transient, yielding occurs 
at a lower yield stress , point (b), corresponding to the elevated 
material temperature. The transient peak strain corresponds to point 
(c) on the curve. 


r 


i 
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Now, as the temperature difference is relaxed and steady state is 
achieved, the stress relaxes along the elastic line from (c) toward (d). 

In fact, for the particular element being examined, the strain relaxation 
is sufficient to produce a very slight amount of tensile yielding at the 
hot conditions (such as the point d'). The cooldown transient results in 
the tensile state, point (e). Cycling continues between these temperature 
limits; and, it is the magnitude of the stable strain range, such as 
between points (f) and (g), which will later be used to predict fatigue life. 

Stress-strain cycles are shown in Figure 11, 12, and 13 for elements 
15, 17, and 30 respectively. Behavior similar to element 34 can be 
identified. The convergence to a stable stress strain cycle is 
demonstrated in Figure 14. The peak strain values for each of the four 
elements 15, 17, 30, and 34 is plotted versus the calculation loop 
number. Note, that a constant value of strain range is approached for 
each element. 

Sample output from the RETSCP program is given in Appendix C. Data 
from the hot and cold transient overshoots is given since these are 
limiting points of the strain cycle. The nodal point displacements 
(inches) and the yield data summary (inches/inch or psi where applicable) 
are given for both parts of the cycle. 



Calculation loop numbers 



Figure li* Effective stress-effective strain diagram for element 15. 





Calculation loop numbers 



Figure 13. Effective stress -effective strain diagram for element 30 





The strain range is determined as follows: the difference between the 
total and plastic strains for the cooling cycle is the elastic strain, 
corresponding to the line 3-4 at the top of Figure 14. Thus, the 
effective strain range for the element, Is the total strain from the heating 
half cycle 4-5-6, plus the cooling elastic strain 3-4. 

The effective strain range distribution, at the combustion chamber 
nozzle throat, is presented in Figure 15. It is obvious that large strains 
in the copper liner are associated with regions of high temperatures. 

Also, as expected, the nickel closeout jacket behaves as an elastic structure. 

Note , that the peak strain range location does not coincide with the 
peak temperature location. That is apparently due to the effect of 
coolant-to-chamber pressure differential. From elastic plate theory, 
pressure loading at run conditions would increase compression on the 
coolant side and relieve compression on the chamber side, along the 
line (i) - (ii) in Figure 15. Thus, the incompressibility of plastic strains 
increases the strain magnitudes toward the liner rib along the hot wall, 
point (iii) in Figure 1 5 . 

By extrapolation (cross plots of the strains at the element centers) to the hot wall 

surface, the maximum strain range is 2.46%, at point (iii) in Figure 15, 

o o 

where the maximum copper temperature is 970 F (521 c). Strain range 
values are generally between 2.0 and 2.5 percent throughout the hot 
portion of the copper liner. 
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Figure 15. Effective strain 



strain range 


distribution (%) at noz- 'e throat 



The strain range results, which are the basis of the life estimate, 
depend on the extremes of the operating cycle. Intermediate points 
of the cycle may, in certain cases, be equally important in determining 
engine life. It was previously noted, that stress relaxation following 
the start transient, could lead to tensile hold strain under hot conditions. 
This phenomena is quite damaging according to the strain range 
partitioning method, Reference 14. Damage of this type is related to 
run time. Reference 19. For this particular case, run times are small 
(total accumulated run time less than 80 seconds) and the tensile hold 
strains are also small. Thus, the important point here is that the 
program has the ability to predict such hot tensile hold strains. 

As an example, distributions of the hoop stress component within 
the copper are shown in Figure 16. The extremes of the cycle 
corresponding to the hot and cold transients are shown by the solid 
lines in Figure 16. Stress distributions at steady run conditions are 
indicated by the shaded area under the dash line curves. As anticipated, 
the stress relaxation is sufficient to produce hot tensile hold states. 
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Chamber outline 



Figure 16. Hoop direction stress component distributions. 



FATIGUE LIFE ANALYSIS 


i 

i 


Prediction of chamber i. • jue life is made by using the computed 
maximum effective strain range value in conjunction with uniaxial 
fatigue test data. A fatigue damage law is used to relate the actual 
operating cycle to uniaxial test conditions. 


Material Properties 

Typical low-cycle fatigue experimental data for OFHC copper, at 
room temperature is shown by the open circle in Figure 17 , Reference 12. 
This data is compared with the well known Universal Slope? life 
predication (c.f. Reference 15); namely. 


4c t = 


3.5 a u _ - 0.12 


n (i-Ra| 


0.60 -0.60 


N 


f 


( 2 ) 


where, the ultimate tensile strength ( au ), modulus of elasticity (E) 
and reduction in area (RA) are based on the short term tensile test. 

The comparison not only shows the applicability of the Universal Slopes 
method; but, also indicates the degree of scatter which is associated 
with fatigue test results . 


Now, for elevated test temperatures, the Universal Slopes prediction 
is modified. Reference 16. First, equation (2) is based on the short 
term properties at the elevated temperature. Then, to account for 
reduced life at elevated temperature, the "average" life is defined as 
one-fifth of the value In equation (2), N^ /5. This procedure was based 
on comparison with a great deal of experimental data. Universal Slopes 
average life is compared, in Figure 17, with isothermal fatigue test data 
for OFHC copper at 1000°F (538°C). 
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Figure 17. Isothermal fatigue of annealed OFHC copper . 
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The above sets of data are based on isothermal fatigue. That is, 
mechanically produced strain at a constant temperature. Thermal 
fatigue test data is obtained by thermal strain cycling between two 
temperature limits. This type of data was obtained for a zirconium- 
copper alloy, which was thermally cycled between 500°F (260°C) 
and 1000°F (538°C) Reference 17. For these conditions, fatigue 
life was lower than for isothermal testing at either of the temperature 
extremes. Other factors may have been involved in this life degradation. 
For the present work , the OFHC copper chamber fatigue life will be 
based on the 1000°F (538°C) data shown in Figure 17 (dash line curve)* 


Fatigue Damage 

The linear cummulative damage law , Reference 2 , is useful when the 
loading produces varying strain range values . According to that law , 
if nj and n^ cycles are applied at conditions for which the fatigue life 
would be and respectively; then, the damage is accumulated 
linearly as follows: 



(3) 


Failure occurs when the cumulated damage totals unity. 


For the case at hand, mean or typical loading conditions are assumed. 
Thus, the fatigue life is obtained directly from the material fatigue 
test data . Had the strain calculation been performed for various loading 
conditions; the life would be computed using equation (3). It was also 
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noted that, for the copper test chamber, experimental run times were 
small. Thus, creep or hold time effects are negligable. When this 
is not the case, a method for incorporating hold time effects, such 
as the strain range partitioning method. Reference 14, is employed. 

Much work has been done in the area of fatigue damage, along both 
the experimental and analytical lines . The topic is touched on 
briefly here only to indicate how an accurate loading spectrum would 
be treated. Also, a precautionary remark is included, that engine 
firing hold time is critical to the operating life. 

Results 

Recall that the maximum value of effective strain range for the 
twentieth test cycle of the copper chamber was 2.46% at the hot 
wall centered under the rib. Based on the isothermal fatigue data 
at 1000°F (538°C), dash line in Figure 17, the predicted fatigue life 
for the copper chamber is 80 cycles. The point represented by the 
copper chamber life (39 cycles), and the maximum strain value for 
the twentieth cycle, is indicated by x in Figure 17. 



CONCLUDING REMARKS 


A detailed life analysis was conducted for a configuration and loading 
conditions which correspond directly to actual thrust chamber test 
conditions. In this manner a direct comparison of predicted and 
observed fatigue life is possible. 

The analysis approximates the operating cycle by a succession of 
quasi-steady states. Thus, actual points within a cycle, and cycle 
to cycle variations can be identified. For example, this is the first 
known computation of tensile hold strain at hot engine test conditions. 

Hold time strain damage analysis was not performed since total hold 
time at steady state run conditions was less than 80 seconds for this 
particular chamber. However, for other configurations with longer 
test firings , hot tensile hold periods could be quite damaging as predicted 
by the strain range partitioning method. 

The mathematical representation of the OFHC copper loading spectrum 
could be improved as follows: Calculate the strain-cycle at thermal 

conditions corresponding to the start and end of the test program; for 
example, cycles five and forty. Account for cycle by cycle strain 
hardening (isotropic or kinematic hardening) until hardening is sufficient 
that cyclic stress-strain data is applicable. Use the linear-cummulative 
damage law with the individual sets of strain data to predict life. 
Geometric changes such as coolant passage bulging, and the bulging 
effect on heat transfer, could also be included in such an analysis. 


Finally, the full three-dimensional capabilities of the RETSCP program 
could be exploited. For example, axial curvature at the throat, with 
the resulting change in boundary constraint, could be modeled. 

Also, axial temperature distribution could be included. It is suspected 
that axial strain relief would be indicated by both of these three- 
dimensional phenomena. 

It is the author's opinion that analytical refinements, such as those 
extended in this report and those proposed in the concluding remarks , 
not only produce more accurate life estimates; but also lead to greater 
understanding of the structural behavior of reusuable rocket engines. 



APPENDIX A — SYMBOLS 


D 

E 

m 

m o 

N f 

P 

P 

c 

RA 

T 

t 


A r 


o 

V 

Ai 

A, 


Fatigue damage fraction 
Modulus of elasticity 
Plastic modulus ratio 

Plastic modulus ratio at reference temperature 
Universal slopes cycles to failure 
Number of cycles to failure 
Pressure 

Chamber pressure 
Reduction in area 
Temperature 

Thrust chamber run time 
Cartesian coordinates 
Total strain range 
Strain (effective) 

Stress (effective) 

Yield stress 

Yield stress at reference temperature 
Poisson's ratio 

Yield stress slope with temperature (Aa^/AT) 
Plastic modulus slope with temperature (Am/AT) 



APPENDIX B — SAMPLE RETSCP INPUT DATA 


The following input listings are from the hot transient overshoot 
portion of the test cycle. The data are input from the sixth 
calculation loop. The data are included as programming guidelines. 
For complete description of input data and format see Reference 1 . 
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APPENDIX C SAMPLE RETSCP OUTPUT DATA 


The following data from the hot and cold transient overshoots is 
given since these are limiting points of the strain cycle. The 
nodal point displacements (inches) and the yield data summary 
(inches/inch or psi where applicable) are given for both parts 
of the cycle. 
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Hot Transient Overshoot - Condition (b) *- Sixth Calculation Loop 
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Cold Transient Overshoot - Condition (d) - Sixth Calculation Loop 
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